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1. Introduction 

One of interesting features inherent in the theories with more than four spacetime dimen- 
sions is a possibility to explain the mysterious pattern of fermion mass hierarchies [||, ||, |3| 
(see also review M and references therein). In the previous series of works 11,1,11, we 
have constructed a model in which a single family of fermions, with vector- like couplings to 
the Standard Model (SM) gauge groups in 6 dimensions, gives rise to three generations of 
chiral Standard Model fermions in 4 dimensions. This mechanism is based on localization 
of fermion zero modes on a two-dimensional vortex with winding number n. In particular, 
in the paper Q the case of global vortex (i.e. group of symmetry of the vortex Ug{l) is 
global) with winding number n = 3 was considered. In this setup three localized fermionic 
zero modes appear due to specific Yukawa coupling to the vortex scalar ^. Coupling of 
fermions to the SM Higgs doublet H results in four-dimensional effective fermion masses. 
Inter-generation mixings occur due to explicit breaking of Ug{l) symmetry. The latter 
point is the main drawback of this model since it does not allow to gauge Ug{l) symmetry. 
In the paper this problem was overcome. The price for unbroken Ug{l), however, is the 
necessity to invoke higher dimension operators in the scalar-fermion interactions^. More 
precisely, to obtain n fermionic generations we consider the vortex with winding number 
1 but fermions coupled to the vortex scalar raised to the nth power. In both cases (global 
and gauged Ug{l) symmetry) the hierarchical pattern of the fermion masses occurs due 
to different profiles of the fermionic wave functions in the transverse extra dimensions. 

^The requirement of renormalizability of the theory does not make much sense anyway in the six- 
dimensional models, since even usual Yukawa scalar-fermion-fermion coupling is non-renormalizable. 
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Namely, in Refs. |§, |3| a crude dimensional analysis has been presented which shows that 
hierarchical mass pattern has the form 

i:g^:g\ (1.1) 

where g is small Yukawa coupling, which yields the localization of fermionic modes. Also, 
m Refs. @, I it has been pointed out that any exact prediction of the fermion masses (as 



well as a check of the Eq.(l.l)) can be obtained only numerically. 

In this paper we return to the models of Refs. 1^, ^ and pay special attention to a 
check of the Eq.(|l.lD. We will see that this equation literally does not hold. Namely, in 



the case of gauged vortex the hierarchical mass pattern takes the form 

At the same time, in the case of global vortex unnatural hierarchical pattern appears. We 
will also evaluate fermion masses numerically. For reasonable values of parameters in the 
case of the gauged vortex our results reproduce SM fermionic mass pattern. Furthermore, 
we will obtain Cabbibo-Kobayashi-Maskawa (CKM) matrix for the same set of parameters, 
and see that six-dimensional theory reproduces the mixing matrix quite well. 

2. The Setup 

In what follows we restrict ourselves to the model of Ref . [Q| . In this section we give a brief 
description of this construction. Our notations coincide with those used in Refs. |2|, |^. 
In particular, six-dimensional coordinates xa are labeled by capital Latin indices A,B = 
0, . . . , 5. Four-dimensional coordinates are labeled by Greek indices /x, = 0, . . . , 3. 
The Minkowski metric is g^B = diag(-|-, —,...,—). We use also chiral representation for 
six-dimensional Dirac F-matrices (see Ref. ^] ) . Besides this we introduce polar coordinates 
r, and 6 in the x^, x^ plane. 

The matter field content of the six-dimensional theory is summarized in Table ||. There 
are three scalar fields. One of these together with the Ug{l) gauge field, forms a vortex, 
while two other scalars, X and H, develop profiles localized on the vortex. The potential 
term of the Lagrangian which gives rise to the non-trivial profiles for the scalar fields has 
the following form 

V, = ^ (|<I>|2 _ v^f + I (|/^|2 _ ^2)2 ^ ;^2|^|2|^|2 ^ ^ (|^|2 _ ^2^2 ^ ^2|^|2|^|2^ (2.I) 

There is also one fermionic generation which consists of five six-dimensional fermions 
Q, D, U, L and E. These fermions have vector-like coupling to SM gauge bosons, (a subtle 
issue of localization of the latters we do not discuss here), and axial couplings with the 
vortex background^ 

Vi=gg<^'Q^^^Q + gu^*^U^-^U + gd<^*^D^^^D+ 

^i_r i-F (2-2) 

gi<^^L^-^L + g,<^*^E^-^E + h.C. 



^For the scalar-fermion interactions, we take the most general operators of the lowest order, consistent 



with gauge invariance. On the other hand, in the scalar self-interaction (2.1) there could be terms of the 



form $X* or $X*\Hf which we do not include. These terms, as well as terms of the form ^^XQ{1 — Tr)Q 
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fields 




profiles 


charges 


representations 








Ug{l) 


TT ^^\ 


err /'o^ 


err /'Q^ 


scalar 




r [r je 
x* (Uj — U, r [oo ) — V 


1 1 


U 


-1 
i 


i 




scalar 




YlfW — ii„ Yir^^ — n 
vv l^UJ — ^X"! ^ \^CXj J — U 


1 1 


U 


1 


1 


scaiai 


±1 


TJ(r\ 

U(()\ — H(nr;') — D 


1 

— 1 


J-1 /9 
i-i/ Z 


9 




fermion 


Q 


3 L zero modes 


axial +3/2 


+1/6 


2 


3 


fermion 


u 


3 R zero modes 


axial —3/2 


+2/3 


1 


3 


fermion 


D 


3 R zero modes 


axial —3/2 


-1/3 


1 


3 


fermion 


L 


3 L zero modes 


axial +3/2 


-1/2 


2 


1 


fermion 


E 


3 R zero modes 


axial —3/2 


-1 


1 


1 



Table 1: Scalars and fermions with their gauge quantum numbers |^ . For convenience, we describe 
here also the profiles of the classical scalar fields and fcrmionic wave functions in extra dimensions. 



(Fy is a six-dimensional analog of four-dimensional matrix 75). 

Note that the scalar background where the field $ has the winding number one, 
$ = F{r)e^^, has exactly the same topological properties as the background <I>i of the 
vortex with winding number three, $1 = Fi{r)e^^^ . As a result, the interactions Eq. ( |2.2| ) 
give rise to three left-handed (right-handed) zero modes for each of the fermions Q, L (U, 
D, E). All of these modes are localized in the core of the vortex^, and could be identified 
with the three generations of the chiral SM fermions. The last point, however, is not quite 
fair since these zero modes are massless from the four-dimensional point of view. To give 
small (compared to the vortex scale ^/Xv) masses to these modes we introduce the following 
interactions of the fermions with the Higgs doublet H and singlet X 

V2 = YdHXQ^—^D + YuHX*Q^—^U + YiHXL^—^E + h.c, (2.3) 

where Ha = ^apH*^^ a, j3 are SC/t^(2) indices. Finally, there is one more gauge invariant 
in the scalar-fermions interaction 

V3 = YaedH^Q^^^^D + y„e„i?$*QiziZ[/ + YieiH^L^^^^E + h.c, (2.4) 



(we denote Yukawa coupling constants in Eq. ( |2.3| ) as Y^^d^u and in Eq. (2.4) as Y^e^, . . ., 
for convenience). This term yields inter-generation mixings. 



in the (2.2), can be forbidden, for example, by the imposition of Z4, discrete symmetry: $ i$, Q 



exp(ir7^)Q, X X. This symmetry is explisitly broken by the interaction (2.3), (2.4), and, therefore the 



undesirable terms will be generated by the radiative corrections. However, the value of these corrections 
is of order ^^^3 eA^ (for the operator \H\^^X* , for instance; A is the cut off), and for our choice of the 

parameters (see Sec. ^) is of order fO~^ ('S')*-^' where M is mass scale, and A — l/KP. So, if we choose 
the cut off to be of order A ^ lOM, then this corrections can be neglected. 
^ These modes were discussed in detail in Ref. H. 
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It is worth noting that to obtain small fermion masses and mixing, it is not necessary to 
take y or e small. Small masses and mixings originate from small overlaps of the fermionic 
wave functions and the scalar profiles. 

Now let us consider this model in more detail. First of all we will present stable non- 
trivial configuration for the scalars and ^^^(l) gauge fields. Then we will find fermion zero 
modes in this background. Finally, we will calculate the fermion masses and the CKM 
matrix. 



3. The vortex background 

As described in the previous section, in our model there are three scalar fields and one 
Ug{l) gauge field which form a non-trivial vortex background. In this section we present a 
stable solution of the static field equations which meets all the requirements. 

The set of the static equations to the background has the following form (prime denotes 
the derivative with respect to r; e is Ug{l) gauge charge), 

F" + -F' - 4tF(1 - Af - Xv'^F (F^ - l) - h^FH^ - rfFX"^ = 0, 

A' - -A' + 2e'^v^F'^(l -A)- 2e^H'^A - le'^X^A = 0, 
r 

H" + -H' - \hA^ - h^Hv'^F^ - kH{H^ - fi^) = 0, 



X" + ^X' - ^XA^ - TfXv^F"^ - pX{X^ - vi^) = 0, (3.1) 



with the boundary conditions 

F{r) = 0, H'{r) = 0, X'{r) = 0, A{r) = 0, at r = 0, 

F{r) 1, H{r) 0, X{r) 0, A{r) 1, at r ^ oo. (3.2) 

Here we use the standard anzatz for Abrikosov-Nielsen-Olesen vortex = 4,5) 

Ai{x) = --^eijXjA{r), <^{x) = vF{r)e'\ 

Ha{x) = 62aH{r), X{x) = X{r). (3.3) 

Let us briefly discuss the choice of the boundary conditions ( ^^ ) and the anzatz ( ^ ) . 
At infinity the fields should tend to their vacuum expectation values (VEV) which are 
determined from the potential (2.1). In general at arbitrary set of the coupling constants 



the potential ( p.ip may admit several minima including those in which X and (or) H are 
non-zero. However, there are several regions in the space of the parameters where the 



potential (2A) has only one minimum ^ = v, H = X = 0. In particular, one of the regions 
is determined by the conditions h"^ > Xk, Xp > r]'^, ifv^ > pvf, and in what follows we will 
assume that these conditions are satisfied. 

In general, any configuration of the fields <I>, H, and X may be characterized by three 
topological charges Qt = (Q*, Qt ^ Qt)- The zero VEVs of X and H admit the existence 
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of the topologically trivial configuration of these fields (|3.3D (in spite of the interaction with 
winded gauge field), Qt = (1,0,0). In principle, there may exist a solution of the field 
equation with Qt = (1, —1, 1) and with the same topological properties of the gauge field 
A^. However, as we will see, the solution with Qt = (1,0,0) is stable, and, therefore, the 
choice of the anzatz ( p. 3D is apropriate. 

As mentioned in the previous section, scalar- fermion interactions ( |2.3D , ( |2.4| ) generate 
low energy mass terms for fermionic modes. Thus, to obtain non-zero fermion masses one 
should have a non-trivial solution for H{r) and X{r). On the other hand, the boundary 
value problem ( |3.1| ), ( |3.2| ) always adopts trivial solutions H = Q oi (and) X = 0. For- 
tunately, it has been pointed out in [P, that in some region of parameters these trivial 
solutions are unstable. Namely, the necessary conditions to have non-trivial stable solu- 
tions are h?v'^ > KpP' and rfv'^ > pvf Q. However, it turns out that these conditions are 
insufficient. To clarify the situation, let us study the problem of the stability in more detail. 
The stability of the scalar <I> and gauge fields is guaranteed by topology. To investigate 
the stability of the static solution H{r) one linearizes equations of motion ( |3.1] ) and obtain 
Schrodinger-type problem with the potential 

h\'^F\r) + ^A^ - k{^? - 3H^{r)). (3.4) 



At large r the solution of Eq.(3T), H{r), has asymptotic H ~ exp(— y^/i^v^ — Kfi'^ r). 
So, at large enough h the Higgs profile becomes narrow compared to F{r). This means, in 
particular, that there is a region of r where H{r) becomes negligible while A(r) and F{r) 



are still almost equal to zero. Therefore, the main contribution to the potential (3^) in this 
region comes from r-independent term — k//^ and is negative. Thus, the potential becomes 
not positively definite, and there may exist negative levels which mean the instability. The 
numerical calculations show that this is the case. 

To be specific, we have chosen the following set of the parameters: v = fi = 1 ■ M^, 
vi = 0.3 • M2, X = k = p = 1 ■ M-2, h = 1.02 • M-\ r/ = 0.31 • M-\ e = 0.04 • M'^ 
where M is a unit of mass which we will define in Section (^). In numerical calculations 
we assume M = 1. To solve the boundary value problem ( |3.1D , ( |3.2| ) we have used the 
relaxation method [^. The solutions for profiles F{r), H{r), X(r) and A{r) in the vicinity 
of the origin are presented in Figure |l[ The potential (3^) for these solutions is positive 
and, therefore, the solutions are stable. 

The vortex background in the case of global Ug{l) group can be obtained from the 
Eqs.(|]l|), (U) by the setting e = 0, A{r) = (A(oo) = 0). The profiles for the scalar 
fields have the same shapes as in the gauge case. We do not present the corresponding plot 
here. 



4. Mass hierarchy: cases of global and gauged Ug{l) 

Now we have all ingredients to find fermionic zero modes in the vortex background and to 
obtain hierarchy of the low energy fermionic masses. To do this let us first briefly remind 
the analysis given in Refs. [||, ^. 
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Figure 1: Profiles of tlie baclfground a.tX — k — p — l,v — fi — 1, e — 0.04, h — 1.02, ij — 0.31, 
vi = 0.3. 



4.1 Global Ug{l) 

Let us consider first the case of global Ug{l) group and one fermionic family, say Q and 
U. The part of the Lagrangian describing the interaction of six-dimensional spinor Q with 
vortex scalar <1> is 



1 



-Q + h.c, 



(4.1) 



Lq = igr^'dAQ - \y9q^-'Q 

This Lagrangian is invariant under Ug{l) rotations 

In the vortex background, one can expand an arbitrary spinor in eigenfunctions of the 
corresponding Dirac operator. There is a set of massive heavy modes (with masses of 
order gv'^) and zero modes. The careful analysis given in |Q |2|, ^] shows that in the 
vortex background with winding number one there are exactly three localized zero modes 
describing four-dimensional left-handed spinors. To obtain three right-handed zero modes 
one should consider a six-dimensional spinor U which has an axial charge —3/2 under Ug{l) 
with the Lagrangian 



Lu = iUT^dAU - ( gu^'^lJ^^^^ 



[/ + h.c, 



(4.2) 



These left-handed and right-handed zero modes have the form 

/ \ /'U4-p(r)e^(i-P)''\ 

gp(r)e^(3-p)^' 

g4-p(r-)e*(i-P)'' 



\ 







Ur, 



J 








V Up{r)e 



(4.3) 
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Here p = 1,2,3 enumerates three modes (three fermionic generations). The radial functions 
Qp are solutions to the following set of the equations. 



(4.4) 



drQi-p - q^-p + 9qV^F^qp = 

with the normalization condition 

oo 

Jrdr{ql+ql_p) = ^. (4.5) 


The functions Up{r) satisfy the same equations with the replacement Qq ^ Qu- 

Since the background radial function F(r) is not known in analytical form, Eqs.([4.4) 

can be solved only numerically. We will present the numerical results hereinafter, now, 

however, let us consider the equations in more details. 

We are interested in obtaining low energy fermion masses which originate from the 

integration over extra dimensions of the following expression 

y„ / cfxHX*Q^—^U + h.c. (4.6) 



Substituting here the zero modes ( |4.3| ), and integrating over polar angle 9 one finds 
m;, = Yu I rdrdeH{r)X{r)qp{gq,r)us{gu,ry'^'-'^ = 



OO 

5ps2TiYu JrdrH{r)X{r)qp{gg,r)up{gu,r). (4.7) 



It is worth noting that due to orthogonality of the zero modes and due to the fact that 
H and X do not depend on the integral over leads to a selection rule, rripg ~ 5ps- To 
obtain non-trivial inter-generation mixings one should also consider the term 

y„e„ j d^xH^*Q^—^U + h.c. (4.8) 

The non-trivial ^-dependence of $ ( p.3[ ) leads to the off-diagonal elements in the mass 
matrix 

J.e{r-s-l) 



= j rdrdeH{r)^{r)qp{gq,r)us{gu,r)i 



CXJ 

(5p,s+i27ry„e„ j rdrH{r)^{r)qp{gq,r)up-i{gu,r). (4. 
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In the same way one finds non-zero off-diagonal elements for down-type quarks m^^ 
(5p+i^s. Therefore, the mass matrices have the following form 

/mil \ / mil mi2 



m 





771-21 m22 
771-32 



777- 



V 









m22 m23 

77133 



To evaluate the integrals (4.7) and (4.9) over r we take into account that the functions 
H{r) and X(r) are localized inside the vortex. So, the integral (4.7) is saturated near the 
origin, and therefore only the behavior of the zero modes at r — > is relevant. It is easy to 
find from equations (4.4) that qp ^ r'^^P at r ^ 0, however, it is not so easy to obtain the 
dependence on Qg of the zero modes. To obtain the dependence on small Yukawa couplings 



g of the fermion masses we should take into account the normalization condition ( [4.5]) and 
investigate the behavior of the zero modes near the origin more carefully. To do this let us 
consider the following simplified model. 

Let us approximate the vortex background by the functions 

v'Firf 

r > 1 



H{r) = X{r) 



(4.10) 



1, r < 1 
0, r > 1 

These functions have the correct behavior at the origin and at infinity. So, this is a 
reasonable approximation. The zero mode equations ( [4. 4)) with the background Eq.( 4.1C| ) 
have analytical solution in terms of modified Bessel functions. However, to clarify the 
picture, we present here only the leading dependence on g of the solutions at r < 1 (which 
is only relevant for our purposes), and exact solutions at r > 1: 



qi[r) 



Cir^ 

Cr , 



\ CrnT 



r < 1 
r > 1 

Cs, r < 1 

^f,-9qr(i I J_) I 



r < 1 
e-9'?'", r > 1 



(4.11) 



To find the dependence on gq of the coefficients C and C let us assume that the 
normalization integral (4^) is saturated at r > 1. Then substituting the solutions ( 4.1l| ) 
into the condition (|4.5D one has C ~ Cm — gq ■ Thus, from the continuity of the solutions 



at r = 1, we find that Ci^2 — gq and C3 — gq. This means in particular that our assumption 
is indeed valid: the normalization integrals at r < 1 are of order g^ and can be neglected. 
Now we can find the mass pattern in this simplified model. Substituting the solutions 



( 4.11 ) and background ( 4.10 ) into integral (4/7) we have the following dependence on gq^^ 
of the fermion masses, 

mil ~ 77722 ~ (c/g^u)^, 77733 ~ (555„). (4.12) 

We see that 77111 and 77722 are parametrically of the same order, and therefore there is no 
hierarchy in the case of global Ug{l). This result was confirmed by the numerical evaluation 
of the zero modes and the integral 
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4.2 Gauged Ug{l) 

Let us concentrate now on the case of the gauged Ug{l). In this case the equations for the 



zero modes can be obtained from ( [4.4[ ) by the replacement 

3A{r) 



dr dr + 



2r 



Therefore, ah dependence of the zero modes on the gauge field can be absorbed into a 
factor 



(gp)gauge = exp j ^^T^dr'j ■ (qp)^ 
Approximating the gauge field as 

A{r) : 



;lobal • 



0, r < 1 

1, r > 1 



(4.13) 



(4.14) 



we find the following zero modes 



qi{r) 



r < 1 



92 (r) 



Car, 



Cm 



r < 1 

g-9?»- 7- > 1 



r < 1 



fg^ y gqrn — 



(4.15) 



The situation changes drastically compared to the global f7g(l) group. The continuity 
at r = 1 results in the following relations between Cp and C, Cm- 



(4.16) 



Using the normalization condition ( [4. 5] ) one obtains that the main contribution to the 
integral ( |4.5| ) for the second mode comes from r > 1, and, so. Cm ~ Qq- However, the 
normalization integral for the first and third modes is saturated now in the core r < 1 by 

3 /2 

the third mode. Thus, C ~ and 



3-p 



(4.17) 



With the dependence ( 4.17| ), one finds the hierarchical mass pattern in the case of gauge 
vortex 



^pp ~ yuigugqf^" 

which was announced in the Introduction. 

In the same way one estimates mixing terms from the Lagrangian (I 



mps ~ 



^p,s+l^u(^u 



3-P 



and 



QuQq 



o\ 

eu^/g^ 1 y 



(4.18) 



(4.19) 



(4.20) 
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Figure 2: The f/-dependence of the fermion mass pattern rhps = nips/g''^ ^ ^^^"^ calculated in the 
background shown at Fig. |l|. 

Of course, the vortex background obtained numerically differs from our crude approx- 
imation. In particular, one can see from Fig. [l| that the gauge field A{r) is wider than the 
scalar F{r), and even far outside the core is almost equal to zero. So, our approximation 
( |4.14 ) is not quite appropriate. One may expect that in the real background (Fig. |l|) the 



g'-dependence of the mass pattern is more resembling with the dependence in the case of 
the global Ug{l) group ( |4.12| ) when A[r) = 0. However, the numerical calculations show 
that the mass pattern depends weakly on the width of the gauge field. We present in Fig. ^ 
five non-zero elements of the mass matrix ( 4.20| ) divided by the corresponding power of g 



in the case g = gq = gu- We see that irips/g^^ ^ ^^^^ indeed depend weakly on g in the 
wide interesting range. 

5. Results for Standard Model fermions 

Now we are ready to find the set of the parameters of the our model (Y, e, g) which 
reproduce the real SM fermionic mass pattern. 

To do this, one has to find eigenvalues of the mass matrix mpg{Yu, eu,9q,gu), to equate 
the result with known fermion masses, and to obtain three equations on four unknown 
variable, Yu, gg, gu- However, even to solve numerically these equations is not an easy 
problem. The point is that to find the solution one has to obtain the mass matrix rripg 
in some range of the parameters, which requires in turn the calculation of the zero modes 
(solving the Eqs. ( [4.4D ) at each value of g. Fortunately, there is a simpler way. First of all, 
note that the mass matrix is proportional to and, therefore, any ratio of the masses 

does not depend on Forming two ratios of the masses one obtains two equations on 

three variables e^, gg, gu- The second simplification is to use the established dependence 
of the mass matrix on g and e (4.18), ( 4.1S| ). After this the obtained equations can be 



easily solved. In fact, in the quark sector the situation is more complicated since besides 
the mass matrix for up-quarks, there are also the mass matrix for down quarks and the 
CKM matrix. Therefore, we have seven equations (four equations for the mass ratios and 
three equations for CKM angles) on five parameters {gg, Qu, gd, ^u, ^d)- In our simulations 
we have used five of seven equations to reproduce correct mass ratios and one element of 
the CKM-matrix {U§^^^). As a result we have obtained 

5g = 0.0013 • M-^ c/„ = 0.0013 • M~^ = 0.03 • A^-^ e„ = 0.9, = 0.1. (5.1) 

Having at hand these parameters, one can find two other independent elements of the 
CKM-matrix, but, in general, we could not expect that they will coincide with the known 
values. However, we have obtained the following CKM-matrix 

/ 0.9780 0.2084 0.0100 " 
0.2086 0.9769 0.0453 | , (5.2) 
0.0004 0.0464 0.9990 



jjCKM 



and see that this matrix coincides quite well with the measured CKM-matrix in the Stan- 
dard Model g 

/ 0.9742 ^ 0.9757 0.219 ^ 0.226 0.002 ^ 0.005 \ 
U^^^ = 0.219 H- 0.225 0.9734 0.9749 0.037 ^ 0.043 . (5.3) 
(^0.004 ^0.014 0.035 ^0.043 0.9990 ^ 0.9993 y 

In the leptonic sector the correct mass ratios are reproduced at 

5z=5e = 0.01, q = 0.75. 

To find Yukawa couplings Y we should first of all understand how the obtained di- 
mensionful parameters (that is all parameters except e) correlate with the absolute energy 
scale. Let us assume that the first non-zero fermion level is of order lOOTeV. It means that 
gv^ ~ lOOTeV. Substituting here g = 0.001 and v = 1 one finds our units: 

1 • M = lO^TeV. 

With this relation, one can convert all parameters measured in our units to the absolute 
mass scale: 

v = ^i = 10^°TeV^ ui = 3 • lO^TeV^ 

A = K = /> = 10-i°TeV-^ /i = 1.02-10-^TeV-\ i] = 2,.l-lQ-^Te\-^ , e = 4- lO'^TeV"^ 

5g = 5« = 1-3 • 10-2«TeV-^ 5d = 3 • 10"27TeV-^ gi = ge = I ■ l^-^'TeW-^ 

Now taking into account the actual values of the fermion masses one finds the values of 
the Yukawa couplings Y: 

Yu = 3.7 • 10"^■^TeV"^ Yd = 2.1 • 10"2°TeV"^ Yi = 1.8 ■ lO^^^TeV"^. 

It is worth noting that we do not obtain any additional hierarchy: all dimensionless 
combinations of the our parameters are of order unity. For instance, gU^ /y/X = 0.4, 
Yu '^ /gU^ = 1-3, etc. 
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6. Conclusions 



In a class of multidimensional models with one vector-like fermionic family, the low-energy 
effective theory describes three chiral families in four dimensions. Hierarchy of fermionic 
masses appears due to different profiles of the fermionic wave functions in extra dimensions. 
In this paper we have investigated numerically one of the models of this kind. Namely, 
we considered the model suggested in ^ and |^. We have shown that hierarchical mass 
pattern indeed appears in the model |^ with the gauged Ug{l) group which is responsible 
for the forming of Abrikosov-Nielsen-Olesen vortex, but the dependence of the mass matrix 
on the parameters is slightly different than it was presented in the Refs. ^. On the 
other hand, in the case of global Ug{l) unnatural mass pattern appears. In the gauge case 
we found the set of the parameters of the model which reproduces well all known fermion 
masses and mixings, without hierarchy among parameters. Moreover, in SM quark sector 
we have reproduced the nine known values (six masses and three mixing angles) by the 
variation of the seven parameters of the model: three couplings gu,d,q, two eu,d, and two 
Yukawa couplings Yu^d- 
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